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Threshold models of global cascades have been extensively used to model real-world collective 
behavior, such as the contagious spread of fads and the adoption of new technologies. A common 
property of those cascade models is that a vanishingly small seed fraction can spread to a finite 
fraction of an infinitely large network through local infections. In social and economic networks, 
however, individuals’ behavior is often influenced not only by what their direct neighbors are doing, 
but also by what the majority of people are doing as a trend. A trend affects individuals’ behavior 
while individuals’ behavior creates a trend. To analyze such a complex interplay between local- and 
global-scale phenomena, I generalize the standard threshold model by introducing a new type of 
node, called global nodes (or trend followers), whose activation probability depends on a global-scale 
trend; specifically the percentage of activated nodes in the population. The model shows that global 
nodes play a role as accelerating cascades once a trend emerges while reducing the probability of a 
trend emerging. Global nodes thus either facilitate or inhibit cascades, suggesting that a moderate 
share of trend followers may maximize the average size of cascades. 

PACS numbers: 64.60.aq, 89.65.-s. 


I. INTRODUCTION 

The mechanisms underlying various sorts of cascades 
on complex networks have been extensively studied over 
the past decade mm- The threshold model of cascades 
formalized by Watts [lj, among others, has been widely 
employed to describe not only information cascades on 
social networks mm, but also default contagion in fi¬ 
nancial networks PIHII! and cascades on multiplex net¬ 
works mm- 

A common assumption implicitly made in the thresh¬ 
old models is that a node’s activation probability depends 
only on local information: the state of neighbors. The 
point of the Watts model is that a vanishingly small seed 
fraction can spread to a finite fraction of an infinitely 
large network through local infections, suggesting that 
individuals’ behavior at the local scale can add up to a 
global-scale outcome. 

In real-world social and economic networks, nodes (in¬ 
dividuals, firms, banks, etc.) normally collect informa¬ 
tion not only from their direct neighbors, but also from 
a wide variety of media sources. Individual nodes can 
easily access global-scale information (poll results, aver¬ 
age prices, average loan rates, etc.) and use it in their 
own decision-making (whom to vote for, whether to raise 
or maintain prices, whether or not to change loan rates, 
etc.). This implies that a global-scale outcome affects in¬ 
dividuals’ behavior while individuals’ behavior adds up 
to a global-scale outcome. Thus, there arises a complex 
interplay between local- and global-scale phenomena. 

To analyze such an interplay in a formal way, I general¬ 
ize the Watts’ threshold model by considering two types 
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of nodes: local nodes (or trendsetters ) and global nodes 
(or trend followers) . The former are activated only by 
local neighbors as in the standard model, while global 
nodes are nodes whose states are dependent on a global- 
scale outcome. In the model, global nodes will be acti¬ 
vated if a large enough fraction of the total number of 
nodes are activated. 

The main contribution of this paper is to show that 
the introduction of global nodes significantly changes the 
mechanism of how a large-size cascade appears. There 
arises a saddle-node bifurcation in the presence of global 
nodes, and thereby a fold catastrophe may occur. This 
is intuitively because once the global nodes begin to be 
activated, those activated nodes further activates the lo¬ 
cal nodes that are otherwise inactive, leading to a drastic 
cascade. We call this type of drastic cascade the trend- 
driven cascade. 

On the other hand, the global nodes may also inhibit 
cascades, because they are never activated until the cas¬ 
cade size exceeds a certain threshold value. Thus, three 
phases may exist in the model: Phase I: no nodes are 
activated (except for seed nodes); Phase II: only local 
nodes may be activated; and Phase III: every node may 
be activated. It is shown that the boundary between 
Phase II and Phase III is well approximated by a root of 
a fixed-point equation, while the extended cascade condi¬ 
tion proposed by Gleeson and Cahalane [2] can still accu¬ 
rately detect the boundary between Phase I and Phases 
II+III. 

The analysis also reveals that the average size of cas¬ 
cades will be uniquely maximized near the boundary be¬ 
tween Phase II and Phase III, where there are a medium 
fraction of global nodes in the population. These results 
suggest that the share of trend followers in the total pop¬ 
ulation may be a key to understanding real-world collec- 
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tive behavior. 


II. A GENERALIZED THRESHOLD MODEL 
WITH TREND FOLLOWERS 


A. Two types of nodes 


We consider an undirected and unweighted network 
that consists of N number of nodes. The degree of each 
node is denoted by k. Each node is initially “inactive”, 
except for the seed nodes that are originally “active” by 
definition. The seed nodes are selected uniformly at ran¬ 
dom with probability po. The response function for the 
local nodes (or trendsetters) is 


F L (m, k) 


1 iff>0 L 

0 otherwise 


(1) 


where m denotes the number of neighbors that are al¬ 
ready active and $l G (po, 1] is the threshold above which 
local nodes are activated. 

The response function for the global nodes (or trend 
followers ) is given by 


f g (p) 


1 if p > 9q 
0 otherwise 


( 2 ) 


where p is the share of activated nodes in the population, 
and 9q £ (p 0 ,1] is the threshold value. The essential dif¬ 
ference between the two types of response functions is 
that F h is a function of m/k , which is the share of acti¬ 
vated nodes in a vanishingly small fraction of nodes (i.e., 
k/N approaches 0 as A goes to infinity), whereas F G 
is a function of the share of activated nodes in a finite 
fraction of nodes. In this paper, we focus on a particular 
case in which the global nodes are activated if a large 
enough fraction of the total population exceeds a thresh¬ 
old. Note that the value of p is global-scale information 
that is not identifiable from the states of local neighbors. 
The accessibility to global-scale information can be jus¬ 
tified by assuming that each node has access to media, 
such as web sites, TV news, newspapers, etc. 

In numerical simulation, a global node is activated at 
step t if the total share of activated nodes at step t — 
1 exceeds 9q. The local and global nodes are selected 
uniformly at random with probabilities 1 — p, and p. 


B. Treelike approximation for the average size of 
cascades 

A treelike approximation for the conditional size of cas¬ 
cades, denoted by poo, is given by a fixed point of the 


following recursion equations Eiia 

Pt+i = (1 — P)Pt+ 1 + pPt- fi, (3) 

oo k 

p \+1 = Po + (1- Po)J2 Pk J2 B t l (Pt)F h {rn,k), (4) 

k =0 m—0 

P?+i = Po + (1 - po)F G {pt), (5) 

where Pk denotes the degree distribution and B^pf) = 
OT(l — Pt) k ~ m - Pt +1 and p G +1 are the fraction of ac¬ 
tive nodes at step t + 1 among local and global nodes, 
respectively. Note that in Eq.|5b, p G + i is determined in¬ 
dependently of the degree distribution. This is because 
global nodes are assumed to ignore local information 
coming directly from their neighbors, pt is obtained as a 
solution to the following recursion equations 

Pt = (1 - p)p\ + PP ?, ( 6 ) 

oo , k — 1 

p\ = Po + (1 - Po) ~ fc )’ 

k—1 m—0 

( 7 ) 



where z is the mean degree. Thus, six variables 
Pt+i , p\+i , p G |_! , Pt, Pt and P G are determined by six equa¬ 
tions ([3])-([8|, taking p 0 = Po = p G = po as given. No¬ 
tice that poo may not be interpreted as the average size 
of cascades because the probability of a global cascade 
occurring is not necessarily equal to 1. One should in¬ 
terpret poo as the size of the cascade, or equivalently as 
the probability of a randomly selected node being acti¬ 
vated, conditional on the event that at least one seed 
node belongs to the extended vulnerable cluster [3]. 

As discussed by Gleeson El f° r the standard case of 
p = 0, poo becomes an accurate approximation for the av¬ 
erage size of cascades if the initial seed fraction p 0 is not 
vanishingly small, where the frequency of global cascades 
will be sufficiently close to 1. If instead po is vanishingly 
small, then we need to multiply poo by the probability 
that a global cascade occurs. This is because there arises 
a non-negligible possibility that all the seed nodes may 
be located outside the extended vulnerable cluster mu, 
in which case the cascade probability will be less than 
1. In the presence of global nodes, the size of the ex¬ 
tended vulnerable cluster becomes smaller than that of 
the conventional one since global nodes are “immune” to 
local infection. The frequency of global cascades can be 
significantly lower than 1 even if p 0 is not vanishingly 
small. 

The average size of the modified extended vulnerable 
cluster S e is given as 

OO 

Se = ^2pk[l - (i - 9oo) fc ], (9) 

k =0 

oo j 

q t = (1 - p) £ -^[1 - (1 - 1, k). (10) 

L ' z 

k— 1 
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Note that a global cascade occurs if and only if at least 
one seed node belongs to the modified extended vulnera¬ 
ble cluster. The (unconditional) average size of cascades, 
denoted by (S', is approximated as 

S«Poo[l-(l-£e) LP0iVJ ]. (11) 

The terms in the square brackets capture the (approxi¬ 
mated) probability that at least one seed node belongs 
to the modified extended vulnerable cluster. 


III. ANALYSIS 


A. Cascade conditions 

The parameter space within which a global cascade 
may occur is called the cascade region , and the condition 
that determines whether or not a global cascade may oc¬ 
cur is called the cascade condition. The first-order cas¬ 
cade condition is a cascade condition with first-order ac¬ 
curacy while the extended cascade condition guarantees 
second-order accuracy [2]. 

Appendix |A| shows the derivation of the first-order and 
extended cascade conditions. Note that although our 
model is summarized as a system of two recursion equa¬ 
tions, Eqs.( |Al| )-( |A2| ), the only equation relevant to the 
cascade conditions is 

OO 7 

Pt = ( 1 - P) [po + (1 - Po) ^2 — 

k =1 z 

k -1 

( 12 ) 

m—0 

In fact, this recursion equation is the same as that ob¬ 
tained in a hypothetical model in which there were no 
global agents and 100p% of local nodes were randomly 
removed. The cascade region that appears in the present 
model coincides with that of the hypothetical model, 
which implies that it is only the local nodes that can ini¬ 
tiate cascades. This irrelevancy of the global nodes stems 
from the fact that the global nodes are totally unrespon¬ 
sive in the initial stage of contagion under the natural 
assumption of lim /)t _ 1 .j.o j f f^F G (p t - 1 ) = 0. The cascade 
condition in the presence of global nodes is the one that 
indicates the parameter space within which a finite frac¬ 
tion of the local nodes can be activated. 

Fig- 0 demonstrates that our analytical approxima¬ 
tion matches reasonably well with simulations on Erdos- 
Renyi random networks. Fig. [l](a) compares the simu¬ 
lated value of S and the cascade conditions. It is known 
that the first-order cascade condition (dashed line) is not 
very accurate for a finite seed size 13 El, but the extended 
cascade condition (solid line) still correctly predicts the 
upper as well as the lower phase transition. If we set 
Po = P = 0 in the first-order cascade condition, then it 
predicts the cascade region shown in the original Watts 
model [IJ. 
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FIG. 1. (Color online) Percentage of global nodes in the pop¬ 
ulation, p, and average size of cascades on Erdos-Renyi net¬ 
works. (a) Simulated average cascade size (color bar) and the 
cascade conditions. Dashed and solid lines denote the first- 
order and the extended cascade conditions, respectively, (b) 
Lines denote the theoretical average size of cascades, S, and 
markers plot simulated values. N = 10 5 , po = 10/A and 
8l = 8g = -18, averaged over 10 3 random perturbations. 


B. Cascade region for trend-driven cascades 


The commonality in the cascade regions between the 
present model and the hypothetical model with random 
node removal is illustrated in Figj2](a) . The figure shows 
that the conditional cascade size, poo, is identical in the 
two cases for p^ < 9q. For p^ > $g, poo in the 
current model becomes larger than that under the hy¬ 
pothetical model because the “removed” nodes will be 
activated. Note that the size of p x in the presence 
of global nodes increases discontinuously at z, where 
z = inf {2 : poo > 0 q}- This indicates that the parameter 
space can be decomposed into three parts: Phases I, II 
and III [Figs{2](b) and (c)]. Phase I is the region in which 
no nodes are activated (except for the seed nodes). Phase 
II is the cascade region within which only local nodes can 
be activated, while every node may be activated in Phase 
III. 


For Erdos-Renyi random graphs, an approximated 
value of z , denoted by 5, is obtained as a solution to the 
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FIG. 2. (Color online) Trend-driven cascades, (a) Lines de¬ 
note poo obtained from Eqs. while markers plot the 

corresponding value obtained in a model without global nodes 
but 100p% of nodes being removed at random, (b) and (c) 
Colors indicate the simulated S and red solid line denotes the 
extended cascade condition. Blue dash-dotted line denotes an 
approximated lower boundary z above which a trend-driven 
cascade occurs. Roman numerals denote the phase of the re¬ 
gion. N = 10 s , po = 10/IV, 0l = -18, 6q = .18 [in panel (a)], 
averaged over 10 3 random perturbations. 


FIG. 3. (Color online) Bifurcation diagram of the roots p 
of fixed-point equations (p, p) — g(p, p) = 0 [eqs. ©-©]- 
(a) Stability diagram that distinguishes the parameter space 
(z, p) in terms of the number of fixed points. The first-order 
line BC highlighted in pale red denotes a set of combination 
(z,p) on which a trend-driven cascade (i.e., a fold catastro¬ 
phe) occurs. Panels (b)-(d) Bifurcation diagram for alterna¬ 
tive values of p. Black solid and dashed lines denote stable 
and unstable fixed points, respectively. Red line is the re¬ 
alized stable root for cascades starting from the initial seed 
fraction po = 10 -4 . #l = #g = -18. 


following root-finding problem obtained from Eq. (12) 


— z k ~ l e~ s 


(1 - p) [po + (1 - Po) ]T 
k -1 

x ]T B^-\e G )F h (m, k)] - 0G = o, (13) 

m =0 

and subject to z £ [z, z ], where z and z are the lower 
and upper boundaries of the extended cascade region, 
respectively. The upper boundary of Phase III is given 
by z. Note that because p^ and p^ can generally differ, z 
obtained by solving (13) is an approximated, rather than 


the correct, lower boundary of Phase III. Nevertheless, 
Figs j2j[b) and (c) demonstrate that the so-defined z nicely 
matches the numerical lower edge. 


C. Bifurcation analysis and tricritical-point scaling 

The presence of discontinuous transitions can also be 
confirmed by the bifurcation diagram illustrated in Fig. [3] 
As pointed out by Gleeson and Cahalane [2], the high-z 
phase transition (i.e., the second phase transition when 
p = 0 and the third phase transition when p > 0) be¬ 
comes discontinuous due to a saddle-node bifurcation. 
Notice that in this model p must be sufficiently small for 


the phase transitions to occur (approximately p < 0.5 
under the baseline parameters). 

Our particular interest is on fold catastrophes that oc¬ 
cur when the parameter combination (z, p) is on the 
bifurcation curve BC highlighted in pale red, which 
is called the first-order line (Point C is called a cusp 
point T5]). Given p = 0.4, for instance [Fig. dr)], the 
first phase transition occurring at z ss 1.9 (from Phase I 
to Phase II) is smooth, but the second transition (from 
Phase II to Phase III) at z = 2.3(~ z) is discontinu¬ 
ous. The introduction of global nodes distinguishes the 
first smooth transition from the second discontinuous one 
however small the value of p(> 0) is. Note that the classi¬ 
cal cascade region is not identical to the region labeled as 
“1 fixed point” in Fig. [3^a) , because there can be three 
fixed points in phase II. Recently, Lee et al. HU also 
show that a saddle-node bifurcation emerges in a multi¬ 
plex network when there is heterogeneity in the response 
functions. 


Fig-i illustrates the boundary of the classical cascade 
region, called the critical line , detected by finding the 
local maxima of the number of iterations (NOI) in the 
recursion equations (Al)-(A2). It demonstrates that as 
z increases, the NOI changes drastically when crossing 
the critical line and the first-order line. Both lines pass 
through the tricritical point T = (z t ,p t ) = (3.14,0.49) 
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FIG. 4. (Color online) Scaling at the tricritic.al point T. NOI 
denotes the number of iterations needed to ensure that the 
successive iterates in the recursion equations |ATJ-(jA2| differ 
by less than 10~ 10 . The critical line (or the cascade region) 
is identified by the local maxima of the NOI. The critical 
line and the first-order line pass through the tricritical point 
T = (3.14,0.49), where the global maximum of the NOI is at¬ 
tained. (Inset) Scaling of the critical line in a new coordinate 
system (/t z ,k m ) centered at a point on the critical line. The 
log-log plot demonstrates that the critical line in the new co¬ 
ordinate system centered at T (black squares) fits a straight 
line with slope 2.0 near the origin. To obtain the new coor¬ 
dinate system, the numerical critical line is smoothed by a 
7th-order polynomial function. #l = #g = -18. 


(subscript t stands for “tricritical”), where phases I, II 
and III terminates. The NOI attains the global maximum 
at T ITT)) . The tricritical point T and the cusp point C in 
Fig. pfFa) generally differ (recall that the cascade region 
is not distinguished by the number of roots p), but in the 
present case they are very close. 

The inset figure shows the scaling of the critical line in 
the scaling fields mmm, where k z is tangent 

to the critical line at T and is perpendicular to it. A 
new coordinate system is given by the transformation 



cos 9 sin 0 
— sin# cos 9 


z — Zt 

V-tH 


(14) 


where 9 = tt + 9, and 9 (< 0) is the angle of n z with the 
z-axis. In the new coordinate system, the critical line 
near the origin is described by 


A /4>t 


(15) 


where <pt is called the crossover tricritical exponent. In 
the baseline case we have 1 /<j> t ~ 2.0. The log-log plot 
demonstrates that at a point slightly deviated from T, 
the scaling does not follow Eq. (151. 



FIG. 5. (Color online) Cumulative fraction of activated nodes 
in a contagion simulation. A trend-driven cascade occurs 
when the cumulative fraction exceeds 9q. The path labeled 
as “random removal” (red crosses) corresponds to the case 
in which p = 0 but 30% of the local nodes are randomly 
removed, z = 2, N = 10 5 , po = 10/A and #l = .18. 


IV. EFFECT OF TREND-FOLLOWING 
BEHAVIOR ON THE AVERAGE SIZE OF 
CASCADES 


A. Two opposite roles of trend followers 


Fig- m illustrates an example path of the cumulative 
fraction of activated nodes. When z G (z, z). the num¬ 
ber of activated nodes in the steady state can be larger 
than that under p = 0. The figure also shows that the 
speed of convergence is decelerated by the introduction 
of global nodes that are never activated in the early stage 
of contagion. 

It follows from the above results that the effect of in¬ 
troducing global nodes is twofold. First, cascades are 
inhibited by the presence of global nodes. The cascade 
region tends to shrink as the share of global nodes p 
increases, and the speed of contagion decelerates accord¬ 
ingly. Global nodes play a role as “immunized patients” 
that prevent infectious diseases from spreading widely. 
Fig! in Appendix [A] illustrates how the cascade region 
shrinks with p. 

Second, the previous figures also demonstrate that 
the size of cascades may increase with p , provided that 
z G ( z,z ). An intuition is that once the total fraction of 
activated nodes exceeds 9&, all the global nodes are acti¬ 
vated regardless of their degree, meaning that a certain 
proportion of nodes outside the largest connected com¬ 
ponent (LCC) can be activated. Recall that the size of 
global cascades is usually constrained by the size of the 
LCC when p = 0 pdSj. 

In total, whether or not the introduction of global 
nodes increases the cascade size depends on the relative 
strength between the negative and positive effects. 
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that attains Phase III. From Eq.(12), an approximated 
upper bound of p*, denoted as p, is given as 

p = max{0, p}, 


where p = 1 - d G [p 0 + (1 - po) ^ —~ 

k=l Z 

k—1 

( 16 ) 


m =0 



taking the parameters other than p as given. The dot¬ 
ted vertical lines in Fig 6|a) plot p. Another interpre¬ 
tation of p is that it indicates the minimum fraction of 
global nodes that is required to eliminate the possibility 
of trend-driven cascades. A further increase of p from p 
will prevent any global cascades. 

Figure |6][b) compares the theoretical and simulated 
values of p*. The theoretical value of p* is computed 
by a greedy search. It shows that the accuracy of theo¬ 
retical p* gets worse as z increases. This is because the 
larger the mean degree, the lesser the benefit of increas¬ 
ing p , which would undermine the accuracy of the greedy 
search in the presence of a finite-size problem. The aver¬ 
age size of cascades is almost unchanged for 2 > 3 since 
most of the nodes belong to the LCC. 


FIG. 6. (Color online) Optimal value of p. (a) Effect of 
varying p on the average size of cascades. Gray line de¬ 
notes the relative size of the cascade window RW, defined 
as RW(p) = (z— z)| M> o/(s — z)| M= o- Vertical dotted lines de¬ 
note p calculated by ( | 16| ). Other lines and markers show the 
theoretical and simulated values of S , respectively. The max¬ 
imum of S is pointed by arrow, (b) Optimal value of p under 
alternative seed sizes. Lines and markers denote theoretical 
and numerical results, respectively. N = 10 s , 9o = 9 q = .18 
and po = 1CU 4 [in panel (a)], averaged over 10 3 random per¬ 
turbations. 


B. Optimal value of p 

1. Maximizing the average cascade size 

Figure |6][a) illustrates how the average size of cas¬ 
cades S varies with p , taking z as given. It shows that 
S initially increases due to the positive effect described 
above, but after passing a critical value, the negative ef¬ 
fect comes to dominate the positive one and S falls dras¬ 
tically (Gray line shows the relative size of the cascade 
window, defined as RW(m) = (z — z)\p>o/{z — z)| At=0 ). 
Thus, there exists a unique maximizer of S , which we 
call the “optimal” value of p, denoted by p*. The case 
of p* > 0 would be of interest because that is the case in 
which the presence of trend followers facilitates cascades. 
A similar discussion on the maximization of information 
cascades is given by Nematzadeh et al. [5j. 

Although an analytical expression for p* is difficult to 
obtain, we know that p* is within the parameter space 


2. Effect of varying the seed size 

Figj6)(b) also states that po and p* have a positive re¬ 
lationship. To understand this, recall that the increment 
of S can be decomposed into two factors from the rela¬ 
tion S = Poo^, where T = [1 — (1 — S' e )L p ° JV J]. A rise 
in S stems from an increase either in p <*, or in T or in 
both. However, for a given p , the size of poo is deter¬ 
mined solely by the size of the LCC independently of 
po if po -C 1. The only factor that the initial seed size 
can affect is therefore T, the (approximated) frequency 
of cascades. 

Fig-0 illustrates the decomposition of an increment of 
S , defined as AS = S\p =c —S | M= o for c G [0,1]. In the fig¬ 
ure, we see that Ap^ = poo\^=c~ Poo|/j=o is independent 
of po , while AT = Tj /i=c — T| M= o (< 0) increases with p 0 
because the larger the initial seed fraction, the higher the 
probability of a global cascade occurring. Suppose that 
p has increased from 0.2 to 0.4. This largely lowers T 
by shrinking the size of the extended vulnerable cluster 
S e , while increasing p^ by activating more nodes located 
outside the LCC. The former effect dominates the latter 
for a small po, but if p 0 is large enough, then the increase 
in poo would compensate for the reduction in T, which 
could make p = 0.4 optimal. 

V. CONCLUSION AND DISCUSSION 

In this work, I investigated an interplay between local- 
and global-scale contagion by generalizing the standard 
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FIG. 7. (Color online) Decomposition of AS. (a) Changes in 
poo and caused by an increase in p. Arrows indicate the 
points above which a rise in p x dominates the reduction in 
4/, leading AS to be positive, (b) Change in the average size 
of cascades. = #g = -18, z = 2.5. 


threshold model. In real-world social and economic net¬ 
works, it is quite common that nodes respond to certain 
information in diverse ways. For example, a retail firm 
may raise prices if some of its wholesalers (i.e., neighbors) 
increase their prices while others may set prices with ref¬ 
erence to the rate of inflation (i.e., a global-scale trend). 
There is thus response heterogeneity among individuals, 
a topic also discussed in Lee et al. m • Interestingly, al¬ 
though Lee et al. introduced response heterogeneity in 
a multiplex version of the threshold model, they show a 
property similar to one suggested in the present work; a 
global cascade may appear abruptly yet slowly. In the 
context of multiplex networks, the current model could 
be expressed as a 2-layer multiplex model, where local 
and global nodes are located in different layers. The 
“global layer” then forms a complete network where ev¬ 
ery node is connected to every other nodes, while the 
degree distribution in the “local layer” is given by p^. 
In that environment, the results of the current model 
will be replicated if it is assumed that fraction p of the 
total nodes are activated only in the global layer and 
1 — p of nodes are activated only in the local layer. How¬ 
ever, the relationship with the model of Lee et al. [15] 
is not straightforward because they consider two alterna¬ 
tive response functions to capture the “AND” and “OR” 
strategies m- 

Global cascade is not a phenomenon that appears in 
every social and economic network on a regular basis, 
suggesting that many real-world networks are stable. Ac¬ 
cording to the simple model presented here, the observed 
stability of actual networks may be attributed to the 
presence of trend followers who are activated only when 
a large enough fraction of the total population are acti¬ 
vated. The present model offers a possible explanation 
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FIG. 8. (Color online) Extended cascade conditions under al¬ 
ternative values of p. Cascade region in the absence of global 
nodes (p = 0) is painted in red. The dark-red solid (p = .2), 
green dash-dotted (p = .4) and blue dotted (p = .6) lines re¬ 
spectively indicate the outer boundary of the corresponding 
cascade region. 9g = -18, po = ICC 4 . 


of why cascades are not very frequent in real-world net¬ 
works as well as why large cascades may appear abruptly. 
I hope this work will contribute to a better understanding 
of the interplay between local and global phenomena that 
could play a crucial role in the mechanism of collective 
behavior. 
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Appendix A: Cascade condition 

To obtain the cascade condition, we reorganize Eqs. © 
- ([8]) as the following two recursion equations that involve 
two variables, p t and pp. 

oo j k —1 

Pt = {1 - p)[p 0 + {1 - Po)J2 — 5Z 1) 

k —1 m =0 

x F l (to, k)] +p[p 0 + (l- p 0 )F G (p t _ 1 )] , (Al) 

oo k 

Pt = (! - P)[P0 + (1 - Po)J2 Pk B m(pt- 1) 

k—0 m—0 

xF h (m,k)] + p[p 0 + (l-p 0 )F G {p t _ 1 )], (A2) 

































The first-order cascade condition is that the Jacobian ma¬ 
trix of function g, defined as (pt,Pt) = g(pt-i, Pt-i), has 
the largest eigenvalue larger than 1 at (0,0) [2]- How¬ 
ever, since lim Pt l |o dpt 1 F G (pt- 1 ) = 0, the Jacobian 
has rank one, meaning that the trace is the only nonzero 
eigenvalue. Therefore, the first term in Eq.( |Al[ ) is the 
only relevant equation for the cascade condition. 

The first-order cascade condition is given as 

U/0lJ , 

(l-p)(l- Po ) J2 —(*-!)>!, (A3) 

tl z 

The irrelevancy of F G near the origin also allows us to 
borrow the extended cascade condition proposed by Glee- 
son and Cahalane |2j. The extended cascade condition is 
given by b 2 — 4ac < 0, where a, b and c are the coefficients 
in the second-order approximation of Eq. © which is of 
the form aq 2 + bq + c = 0. In the presence of parameter 
p > 0, the extended cascade condition leads to (to first 
order in p Q ) 


or Eq. (A3), where 


oo 7 k —1 

E^E* 


k —1 m =0 

= Cq + C\p + C 2 P 2 , 


(. p)F h (m,k ) 


(AS) 


and 

0 = E E( t 7 1 )(m) ( - 1 )‘ + ”^ f ’ l( ”*' fc )- 

k—l-\-l m=0 ' / V / 

(A6) 

Fig|8] illustrates the cascade region indicated by the ex¬ 
tended cascade condition. 


[(1 — 2p 0 )C 2 — 4C2(po + Go — 2Gopo)](l — p) 2 
— 2Gi(l — p 0 )(l — /i) + 1 < 0, (A4) 
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